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We investigate J/ip transverse momentum distribution in a transport approach. While 
the nuclear modification factor Raa{N p ) at RHIC is almost the same as at SPS, the 
averaged transverse momentum square (j>1) and RAAipt) ar e very different at SPS, RHIC 
and LHC and can be used to differentiate from the J/ip production mechanisms in high 
energy nuclear collisions. 

1. Introduction 

From lattice Quantum Chromodynamics (QCD) calculations, there exists a phase tran- 
sition from ordinary hadronic matter to a new state of matter, the so-called Quark Gluon 
Plasma (QGP), at finite temperature. The J/\l/ suppression has long been considered 
as a probe of the new state produced in high energy heavy ion collisions [ lj. The pri- 
mordially produced charmonia via hard nucleon-nucleon (NN) collisions are subject to 
subsequent nuclear absorption in the initial stage and anomalous suppression in the hot 
and dense medium. The normal and anomalous suppression are indeed observed in heavy 
ion collisions at the CERN Super Proton Synchrotron (SPS) [[2] and investigated in many 
theoretical models [ [3] . 

Different from the J/ip production at SPS, there are a remarkable number of charm 
quarks in the QGP phase produced in higher energy nuclear collisions at the BNL Rel- 
ativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC), and 
the regeneration, namely the recombination of those uncorrelated charm quarks offers 
another source for J / ip production [ H] . Obviously, the regeneration will enhance the J / ip 
yield and alter its momentum spectrum. 

From recently observed J / ip production at RHIC [ [5] , the nuclear modification factor 
Raa{N p ) as a function of the number of participant nucleons N p at RHIC is almost the 
same as at SPS, see Fig{TJ The same suppression at SPS and RHIC looks difficult to be 
understood in models with only initial production mechanism, because the temperature at 
RHIC is higher and then the anomalous suppression is predicted to be stronger at RHIC, 
in comparison with SPS. The puzzle of the same suppression was theoretically studied by 
many models. With the idea of sequential suppression [E], if the temperatures at RHIC 
and SPS are both in between the J/ip dissociation temperature and ip' and Xc dissociation 
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temperature, the value Raa ~ 0.6 will not change from SPS to RHIC energy. Considering 
three gluon fusion as the main J / if) production mechanism [ [7] , the cold nuclear matter 
effect which is medium independent can explain the same suppression too. In the frame of 
regeneration, the competition between the initial production and regeneration can explain 
well the same suppression [ ESI [9] . 

The transverse momentum distribution contains more dynamic information on the char- 
monium production and suppression mechanism. The regenerated J/if>s are mainly dis- 
tributed in low p t and central rapidity region [ [9] , but the high p t region is closely related 
to the Cronin effect [ [10] and leakage effect [ [11] for the initially produced J / if)S. While 
the J /if) yield which is a global quantity is not sensitive to the detailed dynamics, the 
transverse momentum distribution at RHIC is very different from SPS, see FigfSJ and 
may be used to differentiate from the production and suppression mechanisms at differ- 
ent energies. In this paper, we investigate the J/if> transverse momentum moments and 
the p t dependence of Raa from SPS to LHC energy. 

2. Transport Model and Numerical Results 

We start with a transport equation [[121 [9] for the distribution function /^(p t , x t) r|b) 
in central rapidity region and transverse phase space (p 4 ,x t ) at time r and fixed impact 
parameter b, 

dfe/dT + • V/* = -acyfy + /3 9 , (1) 

where $ stands for J/ip, if)' and Xc, = pi/v/pf + m% is the transverse velocity, and 
a*(pt, x t , r|b) and /3*(pt, x t , r|b) are the loss and gain terms representing the anomalous 
suppression and regeneration in the hot medium. Considering the gluon dissociation 
process g + ^ — >c + c, a is the momentum integration of the dissociation cross section [ 
[T5] multiplied by thermal gluon distribution, and f3 can be obtained from detailed balance 
between the suppression and regeneration. The distribution /^(pt, x t , r |b) at initial time 
t is determined by the geometrical superposition of NN collisions, including the Cronin 
effect [[10] and nuclear absorption [ 3J. The local temperature, baryon chemical potential 
and collective flow appeared in the thermal gluon and charm quark distribution functions 
are determined by the ideal hydrodynamic equations [fT2]. 

In comparison with the QGP phase, the particle density in the hadronic phase which 
appears in the later period of the system evolution is much lower. To simplify the numer- 
ical calculation, we neglect the hadron contribution to the J /if) production. Solving the 
coupled transport equation for the charmonium motion and the hydrodynamic equations 
for the QGP evolution, one can obtain the J /if) distribution function at the hadronization 
time and then get the final state J /if) yield and transverse momentum distribution. 

We now calculate the nuclear modification factor Raa and averaged transverse momen- 
tum square (pf) as functions of the number of participant nucleons N p for J/ips produced 
in Pb+Pb collisions at SPS and LHC energy and Au+Au collisions at RHIC. All the calcu- 
lations are in mid rapidity. The corresponding parameters for the initial charmonium and 
charm quark distributions in NN collisions and for the hot medium can be found in Ref.[ 
[9|. At RHIC, both the suppression and regeneration in the medium are stronger than at 
SPS, the competition between the suppression and regeneration leads to almost the same 
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Raa at RHIC and SPS, especially for semi-central and central collisions with N p > 150, 
see FigJTJ However, the case is very different at LHC. The initially produced J/ips are 
almost all eaten up by the very hot, long lived and large fireball, and the regeneration 
becomes dominant in a wide region of N p . Only for peripheral collisions with N p < 50, 
the initial production is important. Due to the increasing suppression and regeneration 
with centrality, the Raa decreases with N p in the initial production dominant region and 
increases with N p in the regeneration dominant region, see FigJT] 




Figure 1. The nuclear modification factor Raa as a function of N p for Pb+Pb collisions 
at SPS and LHC and Au+Au collisions at RHIC at mid rapidity. The SPS and RHIC 
data are from [ [2] and [ [5] , and the lines are our theoretical calculations. 

Now we come to the p t distribution. FigfS] shows the averaged transverse momen- 
tum square normalized by the corresponding value in NN collisions, {pDaa/ {Pt)pp-> as a 
function of centrality (left panel) and the Raa as a function of p t (right panel). While 
the Raa{N p ) is almost the same at SPS and RHIC, the transverse momentum distribu- 
tion is really sensitive to the production mechanism. The dominant initial production, 
Cronin effect and leakage effect at SPS lead to an increasing (pf) aa/ (Pt)pp with N p and 
an increasing RAAiPt) with p t . In contrast to SPS, the regeneration is the dominant pro- 
duction mechanism at LHC which results in decreasing (pf) aa/ '(Pt)pp anc ^ ^AAipt)- Since 
we assumed charm quark thermalization in the medium, (pf) aa/ ipDpp becomes saturated 
when the contribution from the initial production disappears. Considering the fact that 
the regenerated J ftps carry low momentum, the Raa{pi) is larger than 1 at low pt but 
vanishes at high p t . The case at RHIC is in between SPS and LHC where both the initial 
production and regeneration are important and the competition between them controls 
the J /if) production. It is the contribution from the regeneration that separates clearly 
the J /if) transverse momentum distribution at RHIC from that at SPS. 

In summary, we studied the contribution from the primordial production in the initial 
stage and the regeneration in the medium to the J /if) production. We found that with 
increasing fraction of the regeneration from SPS to LHC, the transverse momentum distri- 
bution behaves very differently and can be used to differentiate from the J / if) production 
mechanism in high energy nuclear collisions. 
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Figure 2. The normalized averaged transverse momentum square (pI)aa/ (Pt)pp as a func- 
tion of N p (left panel) and Raa at 6 = as a function of p t (right panel) for Pb+Pb 
collisions at SPS and LHC and Au+Au collisions at RHIC at mid rapidity. The SPS and 
RHIC data are from [[H] and [US], and the lines are our theoretical calculations. 
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